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I. INTRODUCTION 

In a recent paper]] hereafter denoted by I, a local 
field theory capable of describing the zero temperature 
(T = 0) ferromagnetic phase transition in disordered itin- 
erant electron systems was developed. In the present pa- 
per this theory is used to exactly determine the critical 
behavior at the phase transition, and the connections be- 
tween the local theory and previous descriptions of the 
ferromagnetic quantum phase transition are established. 

Historically, the ferromagnetic transition in itinerant 
electron systems at T = was the first quantum phase 
transition to be studied in detail. Hertzcl concluded 
that the transition in the physically interesting dimen- 
sion d = 3 was mean field-like. The basic idea behind 
this result was that the effective dimension of the sys- 
tem, which is given by the spatial dimension d plus the 
effective time dimension z, was above the upper criti- 
cal dimension for the transition, so that fluctuation ef- 
fects could be ignored. This conclusion is now knpwn 
to be incorrect. For example, the Harris criterion!] for 
phase transitions in disordered systems states that the 
correlation length exponent v must satisfy the inequal- 
ity v > 2/d, whereas mean- field theory gives v = 1/2 
for all d. This in turn implies that a simple mean-field 



description must break down in dimensions d < 4. 

The reason for this breakdown of mean-field theory 
was shown in Ref. ^ to be the existence, in itinerant elec- 
tron systems, of soft or massless modes other than the 
order parameter fluctuations, which were not taken into 
account in Hertz's theory. In disordered systems these 
modes are diffusive, and they couple to the order pa* 
rameter fluctuations and modify the critical behavior 
Among other things, they lead to a correlation length 
exponent that satisfies the Harris criterion. Technically, 
if these additional soft modes are integrated out, they 
lead to a long-ranged interaction between the order pa- 
rameter fluctuations. It was argued that once this effect 
is taken into account, all other fluctuation effects are sup- 
pressed by the long-range nature of the interactions and 
that the critical behavior is governed by a fixed point that 
is Gaussian, but does not yield mean-field exponents. It 
was thought that the critical behavior found in Ref. || 
was exact. 

Several years before the work reported in Ref. |J, the 
study of metal-insulator transitions of disordered inter- 
acting electrons constituted a separate development in 
the many-electron problem.Q Within this context, a tran- 
sition was encountered that was not a metal-insulator 
transition, but rather of magnetic natureJa Due to the 
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methods used in Ref. |], the order parameter and the 
nature of the ordered state were not identified, but the 
critical behavior was determined and found to consist of 
power laws with simple exponents, modified by compli- 
cated logarithmic corrections. The critical behavior for 
the ferromagnetic transition determined in Ref. § turned 
out to consist of the same simple power laws, albeit 
with different, and much simpler, logarithmic corrections. 
This led, in Ref. to the suggestion that the transition 
studied in Ref. |] was in fact the ferromagnetic transition. 
The discrepancy in the logarithmic corrections between 
the two approaches was attributed to the fact that of the 
two integral equations derived in Ref. |8|, only one had 
been shown to be exact. The conclusion thus was that 
the theories presented in Refs. |i| and || had treated the 
same problem, and that the former solution was exact, 
while the latter was approximate. 

The latter conclusion, however, relied on a weak link 
in the chain of arguments, since the theory developed 
in Ref. |] was not very suitable for determining logarith- 
mic corrections to power laws. The reason was that the 
additional soft modes were integrated out to obtain a de- 
scription solely in terms of order parameter fluctuations. 
The resulting field theory was thus nonlocal, which makes 
explicit calculations cumbersome. Consequently, most of 
the arguments used in Ref. ^ to determine the critical be- 
havior were simple power counting techniques that were 
not sensitive to logarithmic corrections to power laws. 

It is the purpose of the present paper, in conjunction 
with the preceding paper I, to settle the remaining ques- 
tions regarding the relation between Refs. ^ and g, and 
the exact critical behavior, including logarithmic correc- 
tions to scaling, at the quantum ferromagnetic transi- 
tion of disordered itinerant electrons. By using the local 
field-theoretic description of I that explicitly keeps all 
soft modes, we show that Ref. || missed marginal opera- 
tors that lead to logarithmic corrections to the Gaussian 
critical behavior discussed there. Moreover, taking these 
marginal operators into account leads to integral equa- 
tions for the relevant vertex functions that are identical 
to the ones derived in Ref. |[ The current formulation 
can further be mapped onto the one of Ref. ^, which 
shows that the transition found in the latter paper was 
really the ferromagnetic one, and that the results origi- 
nally derived in that reference are exact. 

This paper is organized as follows. In Sec. |TJ we first 
recall the results of I. We then use diagrammatic tech- 
niques to derive exact integral equations for the two-point 
vertex functions that appear in the theory. We conclude 
this section by quoting a previous solution to the se equa- 
tions that is valid at the critical point. In Sec. Ill we show 
how some physical observables in the paramagnetic phase 
are related to these vertex functions. We then develop a 
scaling theory to determine the critical behavior of other 
observables of interest, as well as the critical behavior in 
the ferromagnetic phase. In Sec. [V we discuss general 
theoretical aspects of this paper as well as experimental 
consequences of our results. Various technical issues are 



relegated to several appendices. 



II. EFFECTIVE FIELD THEORY, AND ITS 
SOLUTION 

A. Effective action 

In I it was shown that the effective long-wavelength and 
low-frequency field theory that contains the critical fixed 
point and describes the exact quantum critical behavior 
of disordered itinerant ferromagnets is given by the action 
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(2.1) 

Here M"(k), with components l M"(k), is the fluctuating 
magnetization at wavenumber k and bosonic Matsubara 
frequency fl n = 2itTn, with a a replica label, and the 
field b is defined in terms of the magnetization, 

7-^12 ( k ) — buia 2 (~ Y^ 2 ^n,ni— n 2 ( k ) 
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The labels 1, 2, etc., comprise both frequency and replica 
indices, 1 = (ni,a±), etc. The two-point M-vertex is 
given by 

u 2 (k) =< + ad-2|k| rf - 2 + a 2 k 2 . (2.3a) 

Here the nonanalytic term, proportional to |k| rf ~ 2 , re- 
flects the nonanalytic wavenumber dependence of the 
electron spin susceptibility in a disordered itinerant elec- 
tron system, as has been explained in I. For weak disor- 
der, characterized by a mean-free path I » 1/fcp, with 
hp the Fermi wavenumber, the prefactor ad-2 is of order 
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l/kp£, while a 2 = 0(1). For physical values of the spa- 
tial dimension d, and for asymptotically small wavenum- 
bers, the nonanalytic term dominates the analytic k 2 - 
tcrm. However, for completeness and later reference we 
include the latter, which had been dropped from the fi- 
nal effective action in I. to is the bare distance from the 
ferromagnetic critical point. 

The fermionic degrees of freedom are represented by 
the field q; electron number, spin, and energy density 
fluctuations can all be expressed in terms of the J.?"™- 
These are the additional slow modes mentioned in Sec. 
| above. The frequency labels n > 0, m < of the q 
denote fermionic Matsubara frequency indices, i is a spin 
label (i = and i — 1,2,3 correspond to spin-singlet 
and spin-triplet fluctuations, respectively), and the label 
r (r = 0, 3) serves to write the complex-valued q fields as 
two-component real- number valued fields .□ The fermionic 
part of the action is characterized by the two-point vertex 
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The parameter G = 8/nao is a measure of the disorder, 
with (To the bare conductivity. H — irNp/A, with Np the 
bare single-particle or tunneling density of states per spin 
at the Fermi surface, is a bare quasi-particle density of 
states that also determines the specific heat coefficient. 
Finally, c± , c 2 are coupling constants whose bare values 
are related and given by 
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with K t the spin-triplet interaction amplitude of the elec- 
trons. The replicated partition function is given in terms 
of the action by 



determination of the critical behavior of the field theory 
is reduced to the solution of two coupled integral equa- 
tions that were first derived by different methods in Ref. 



1. Gaussian propagators 

In order to set up a loop expansion we will need the 
basic two-point propagators for the above theory. They 
are determined by the Gaussian action, 
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The quadratic form defined by this Gaussian action has 
been inverted in I. For the order parameter correlations 
we find 
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in terms of the paramagnon propagator, 
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B. Perturbation theory to all orders 

We will now show that the effective action given in 
the preceding subsection can be solved perturbatively ex- 
actly, i.e. it is possible to resum perturbation theory to 
all orders. The basic idea is to first show that the n- 
point vertices for n > 3 are either not renormalized, or 
their renormalization is simply related to that of the two- 
point vertex functions. This in turn implies that exact 
self-consistent equations for the two-point vertex func- 
tions can be derived. The net result will be that the 



The dynamical piece of the paramagnon propagator M. , 
whose structure is characteristic of disordered itinerant 
ferromagnets, has been produced by the coupling be- 
tween the order parameter field and the fermionic degrees 
of freedom. 

For the fermionic propagators one finds 
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Here 2? is the basic diffusion propagator or diffuson. fn 
the limit of small frequencies and wavenumbers it reads 
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Physically, T> describes heat diffusionff 

Finally, due to the coupling between M and q there is 
a mixed propagator, 
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2. 3-point and 4-point vertices 

We now determine loop corrections to the tree-level 
theory. We begin by considering the three-point vertex 
whose coupling constant is denoted by c 2 in Eq. (2.1). 
Diagrammatically the bare three-point vertex is given in 
Fig. |IJ Now consider the one-loop renormalizations of 




FIG. 1. Diagrammatic representation of the bare 
b g 2 -vertex. Dashed lines denote M or 6-fields, and solid lines 
denote q-fields. 

this vertex, which are shown in Fig. @. For scaling pur- 





FIG. 2. One-loop corrections to the vertex shown in Fig. [lj 
Solid lines denote {qq) propagators, dashed lines (66) propa- 
gators, and dashed-dotted lines (bq) propagators. 

poses, we can use simple estimates for the propagators. 
Specifically, the 9-pranagators all scale like an inverse 
wavenumber squared,til 
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The M-propagator at criticality may scale like a num- 
ber, or like an inverse wavenumber to the power d — 2. 
This depends on the scaling behavior of the frequency, 
which can be different in different contexts, as has been 



explained in I and can be seen from Eqs. (2.6c) and ([2.9]), 
respectively: f2 can scale either like |k| d , as in the para- 
magnon propagator, or like k 2 , as in the diffuson. The 
two possibilities therefore are, 
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Similarly, the mixed propagator, Eq. (2.10), scales like 
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If we use an infrared wavenumber cutoff A, we see that 
the integrals that correspond to the diagrams shown in 
Fig. H all scale like A d ~ 2 . That is, the one- loop renormal- 
ization of C2 at zero external wavenumber and frequency 
is a finite number for all d > 2. More generally, an n-loop 
skeleton diagram has n independent wavenumber and fre- 
quency integrals. Diagrams that contain only solid and 
dashed lines contain In (qq) propagators, and up to n 
(MM) propagators. Similar considerations hold for dia- 
grams that contain mixed propagators. The net result is 
that any n-loop skeleton diagram scales like A n ( d ~ 2 \ All 
of these contributions thus amount to finite corrections 
to the bare value of c 2 . By induction it follows that inser- 
tions do not produce singular contributions either. We 
conclude that there are no singular renormalizations of 
the three-point vertex function in the field theory defined 
by Eq. fg~T]:. 

In addition to the renormalization of c 2 , a new three- 
point vertex with a replica structure that is different from 
the one of c 2 is generated by the renormalization group 
at one-loop order. As shown in Appendix the fre- 
quency structure of this vertex is such that it carries one 
more frequency sum and associated temperature factor 
than the vertex with coupling constant c 2 . Its coupling 
constant is therefore more irrelevant than c 2 and can be 
neglected. 

Next we consider the four-point vertex, T^ 4 ) in Eq. 
fl2.ip . The bare four-point vertex is given analytically 
in Eq. ( [2.3c|) and shown diagrammatically in Fig. |j| and 




FIG. 3. The bare four-point vertex 
the one-loop renormalizations are shown in Fig. 0. Notice 





FIG. 4. One-loop corrections to the vertex shown in Fig. |j| 
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that this vertex is proportional either to a wavenumber 
squar ed, o r to a frequency. Using the estimates given by 
Eqs. (2.11), we find that the renormalization of the part 
of r*- 4 - 1 that is proportional to a wavenumber squared is 
always a finite number, again scaling as A^ -2 ^. For the 
part that is proportional to frequency, on the other hand, 
power counting shows that this term can have logarithmi- 
cally singular renormalizations. An explicit calculation 
would be very cumbersome. However, it is not neces- 
sary since the one-loop renormalization of the coupling 
constant H in Eq. (2.3c) obtained this way is identical 
to that obtained by rcnormalizing the two-point vertex, 
Eq. ( 2.3b] ). This is because both terms arise from the 
same term in the underlying nonlinear sigma model for 
the fermionic degcees of freedom, which is believed to 
be renormalizable.113 By the same argument, the renor- 
malization of r^ 4 -* to all orders is given by that of T {2 \ 
and therefore need not be considered separately. The 
explicit calculation of confirms the existence of the 
logarithms that were alluded to above, as we will demon- 
strate in the next subsection. 

In addition to the diagrams that renormalize r*- 4 ), 
there are one-loop terms that represent four-point ver- 
tices with more restrictive replica structures. These cor- 
respond either to the two-body interaction terms that 
were shown in I to not change the critical behavior, or to 
many-body interactions that are shown in Appendix [A] to 
be more irrelevant than I^ 4 ) and thus can be neglected. 



3. 2-point vertices 

We now turn to the two-point vertices in the effective 
action. The one-loop renormalization of the 65-vertex 
is shown in Fig. pi Using Eqs. (2.11), it is easy to see 




+ 



FIG. 6. One-loop renormalizations of the fermionic 
two-point vertex . 



H{in n ) = H + § GK t ^ h £ W Mi(p) 



n n ^ v 

1=0 p 



(2.12) 



which diverges logarithmically as f2„ — > for all 2 < d < 
4. As was explained in I, this divergent renormalization, 
which arises from the nominally irrelevant vertices r^ 4 -* 
and C2 , is a consequence of the presence of two time scales 
in the problem. In addition, there are terms that are fi- 
nite in d > 2. It is importan t to note that the structure 
DM. in the integrand of Eq. (2.12) stems from the second 
term, proportional to T> 2 A4, in the triplet qq-propagator, 
Eq. (2.8b), times a term 2? _1 that is due to the wavenum- 
ber a nd frequency dependence of the quartic vertex, Eq. 

O- 

The exact vertex can be generated from the one- 
loop diagrams by dressing all propagators and all vertices 
in Fig. pi The relevant vertices are and c 2 . As was 



shown in Sec. II B 2, the latter has no singular renormal- 



izations in d > 2, so it need not be dressed. Denoting the 
exact four-g vertex by a square, and the dressed propaga- 
tors by double lines, we therefore have the diagrammatic 
representation of the renormalization of T*- 2 ) to all orders 
shown in Fig. [?]. Analytically, this result corresponds to 





FIG. 5. One-loop renormalizations of the mixed two-point 
vertex. 



that this diagram is finite in d > 2. Since the three- 
point vertex is not singularly renormalized, see the pre- 
vious subsection, it follows that the 6g-vertex has only 
finite renormalizations to all orders in perturbation the- 
ory. This means that the coupling constant K t is not 
singularly renormalized. 

The one-loop renormalizations of r*- 2 ) are shown in Fig. 
As was shown in I, the renormalization of G obtained 
from these diagrams is finite in d > 2. For the renor- 
malized H, which we denote by H(ifl n ), one obtains to 
one-loop order, 



i 



FIG. 7. Renormalization of the fermionic two-point vertex 
r' 2 ' to all orders. 



simply dressing the propagators in Eq. ( 2. 12j ) . Notice 
that this procedure includes the vertex renormalization 
due to the structure pointed out above. Also notice that 
it is crucial for our argument that G and K t carry finite 
renormalizations only. We thus have the exact result, as 
far as the asymptotic critical behavior is concerned, 



1=0 P 



GK t n t + p 2 u 2 (p,in l ) 



(2.13a) 



with u 2 (p,?f2/) the fully renormalized bb- vertex. By the 
same arguments, we obtain the latter as shown in Fig. 
||. By dressing the propagators in the analytic one-loop 
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D. Solution of the integral equations 




FIG. 8. Renormalization of the two-point magnetization 
vertex ui to all orders. 



expression given in Eq. (3.5b) of I, we have, for fl n > 0, 
G 2 Kt ^„ „A 1 



w 2 (k,if2„) = t 



2 V 



1=0 



1 



(p + k) 2 + GH(ini + ffl„)(Qi + On) 



p 2 +GH(iQi)Qi 
(2.13b) 



In writing Eq. (2.13b), we have for simplicity put the 
bare coupling constants aa-2 and a 2 equal to zero, since 
they are generated at one-loop order. 



C. Integral equations for diffusion coefficients 



Eqs. ( 2.13| ) in Sec. II B constitute two closed integral 
equations for the two-point vertices. As we have seen, 
this has been possible to achieve since (1) the four-point 
vertex I^ 4 ) renormalizes like the two-point vertex T'- 2 -', 
and (2) all other vertices are subject to finite renormal- 
izations only. As a result, the solution of Eqs. (2.13) 
provides us with the perturbatively exact critical behav- 



To make contact with previous work, it is useful to 
rewrite Eqs. ( 2.13| ) in terms of the (thermal) diffusion 
coefficient D(in n ) = l/GH(ifl n ) and the spin diffusion 
coefficient, D s (k, ift n ) = u 2 (k, iil n )/Gif t . If we analyt- 
ically continue to real frequencies, i£l n — ► Q + iO, the 
self-consistent one-loop equations read 



D s (k, n) = D° s + — r Y r duJ— 



1 



(p + k) 2 - i(u + tt)/D(uj + n) 



(2.14a) 



1 

DO 



3G 
8V 



y- 



du> 



1 



-iu> + p 2 D s (p,u) 

(2.14b) 



Here D° = l/GH and D° s = t /GK t are the bare dif- 
fusion coefficients. These integral equations were first 
derived in Ref. || by means of a resummation of pertur- 
bation theory within a nonlinear sigma model for inter- 
acting electrons in the limit of a large spin-triplet inter- 
action amplitude. As has been discussed in I, this sigma 
model is recovered from the current model by integrating 
out the magnetization. This is the mapping between the 
two models that was referred to in the Introduction. 



In Ref. |8|, the coupled integral equations, Eqs. (2.14), 
were solved by three distinct methods: A direct ana- 
lytic solution, a renormalization group solution, and a 
numerical solution. The results of all three approaches 
were consistent with one another. Here we will quote 
the most relevant results, restricting ourselves, as in I, to 
2 < d < 4. 

Simple scaling arguments show that at criticality, t = 
M 2 (k = 0,Q„ = 0) = 0, D(il — > 0) is a constant except 
for logarithmic terms. This suggests the ansatz 



D(Q) = D°/F[ln(l/ftT)]. 



(2.15) 



Here 1/r = miG/8m, with n the electron density and m 
the electron effective mass, is the elastic-scattering rate. 
The same arguments yield D s (k,fl = 0) ~ |k| d ~ 2 at the 
critical point except for logarithmic terms. So we write, 

D s (k,n = 0) = D° s (\M/k F ) d - 2 F s [\n(k F /\k\)} . (2.16) 

Solving the resulting equations for F and F s givesi 

D{il -> 0) = D° [ g (lnl/fir))] _1 , (2.17a) 

D s (k ^0,Q = 0) = D° s (\k\/k F ) d - 2 d' (GK t /H) k d F 2 
x [. 9 (21n(fc F /|k|))]- 1 . (2.17b) 



Here 



1 ( X ) = j2i(c(d) X r/i 



,l (n 2 -n)ln(2/d)/2 



(2.18a) 



n=0 



In an asymptotic expansion for large x, the leading term 
is 

3 (x)«[21n(d/2)/7r]- 1/2 e [ ln ( c ^^ 2 / 2 W . (2.18b) 
The dimensionality dependent coefficient c(d) is given by 



c{d)=c'{d)/d'{d) , 



(2.19a) 



where 



d'(d)=c"r(2-d/2)F(2-d/2,l/2;3/2;l) , (2.19b) 

with F a hypergeometric function, T the gamma func- 
tion, and c' and c" smoothly varying functions of d. 

The leading dependence of either D or D s on t away 
from criticality at zero frequency and wavenumber, or 
that of D s at t = as a function of frequency at k = 0, 
will follow from the scaling theory to be developed in Sec. 
[II below. 
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III. CRITICAL BEHAVIOR OF OBSERVABLES 

In this section we determine the exact critical behavior 
of various observables. We do so by developing a general 
scaling description for the free energy and various trans- 
port coefficients and relaxation rates, and using the exact 
solution of Sec. || to determine the values of the indepen- 
dent critical exponents. Again, we restrict ourselves to 
the dimensionality range 2 < d < 4. The results for 
d > 4 in Rcf. || were exact, and the behavior in d = 4 
can be obtained by combining the solution of the integral 
equations from Ref. || for that case with the arguments 
given below. We have also checked some of our results by 
means of explicit perturbation theory, see Appendix 
Parts of the results presented here have been previously 
published in Ref. [yj! 



A. Identification of observables 

We first discuss how to relate the behavior of some 
physical observables of interest near the quantum criti- 
cal point to the solution of the effective field theory given 
in Sec. ||. We consider the electrical conductivity <r, the 
specific heat coefficient 7c, the tunneling density of states 
N, the spin susceptibility Xs, the heat and spin diffusion 
coefficients D and D Sl respectively, the phase relaxation 
rate t~^, and various quasiparticle properties, in partic- 
ular the quasiparticle decay rate Tgp. The magnetization 
m will be obtained from the free energy in Sec. Ill B be- 
low. 

The conductivity, a = 8/nG = l/p, with p the re- 



sistivity, is proportional to the inv erse of the renormal- 
ized disorder parameter, G, in Eq. (2J), and the specific 
heat coefficient, 7c = C/T, with C the specific heat, 
i s proport ional to the renormalized value of H in Eqs. 
(2.3b,2.3c).Q The single particle density of states (DOS) 
as a function of the distance in energy or frequency space 
from the Fermi surface, N(e), is given by Eq. (2.29d) in 
I. In terms of expectation values of the field q, it takes 
the form of an expansion 



N(e F +e) =N f 



(3 m i,r 



, +l0 + O((q 4 }) 



(3.1) 



Experimentally, e is equal to the electron charge times 
the bias voltage. The dynamical spin susceptibility is 
given by the M • M correlation function, and we have 
seen that the coupling between the M and q fields affects 
the dynamical part of that correlation function only. This 
implies that the static spin susceptibility, Xs(k), is given 
by the renormalized value of «2(k, if2„ = 0), cf. Eq. 
( 2.131 j ) . Alternatively, it is proportional to the inverse of 
the static spin diffusion coefficient, 



X.(k) 



1 



2GK t D s (k,m n = 0) 



(3.2) 



The heat and spin diffusion coefficients D and D s are 
explicitly given by the solution of the integral equations 
discussed in Sec. ||. In Appendix ^ we show how to 
generalize the integral equations to the ordered phase, so 
that the equation of stat e or the magnetization, m, can be 
obtained. In Sec. Ill B| we determine the magnetization 
by means of a scaling theory. 

The relaxation rates are defined in terms of the dif- 



fuson, whose bare propagator is given by Eq. (2J)). Its 
renormalized counterpart has the form 

T>{k,iQ n ) 



k 2 + GH{in n )n n k 2 + n n /D(in n ) ' 

(3.3a) 

Here Z is the wavefunction renormalization, which deter- 
mines the single-particle DOS, Eq. (fO|), via N = N-pZ. 
In addition, the quasiparticle DOS Nqp is related to N 



N QP = N/a QP 



(3.3b) 



with oqp the quasiparticle weight. Upon an analytic 
continuation to real frequencies, H(iQ, n ) acquires a real 
part H' and an imaginary part H" , H(ifl n — ► e + iO) = 
H'(e) + iH"(e). Dividing by GH', the renormalized dif- 
fuson can be written 



D(k, iVt r , 



iO) 



(7V QP /iV 2 )a 2 



1 



G £ QP k 2 



QP 



(3.3c) 



'QP 



with a QP = N F Z/H'(e), N QP = H'(e), D QP = l/GH'(e) 
the quasiparticle diffusion coefficient, and 



r QP = eH"(e)/H'(e) 



(3.4a) 



the quasiparticle decay rate. In contrast, the phase 
breaking rate is the "mass" that is acquired by the 
diffus on at real frequencies and finite temperature. Equa- 
tion ( [3. 3a ) shows that it can be identified, apart from a 
multiplicative constant, with 



= Re{SlH{in))/N F 

iU— >e-H0 

or, in terms of the imaginary part of H, 



' P h 



eH"(e)/N F 



(3.4b) 



(3.4c) 



For later reference we also give the renormalized para- 
magnon propagator. Leaving out terms that are irrele- 
vant for our purposes, it reads 



M(k,in r , 



1 



W2 (k,^„) + GA' t |O ri |/k 2 
l/GK t 

D s (k,tn n ) + \n n \/k 2 ■ 



(3.5) 



7 



The above quantities have all been defined at T = 0, 
with an eye on the fact that the results of Sec. || are 
only valid at T = 0. In general, we are also interested 
in the analogous results at T > 0. This is of particular 
importance for the relaxation rates, since only at T > 
do they provide a mass for the diffuson propagator and 
hence constitute a true real-time decay rate. In order to 
obtain complete results at finite temperatures, a Matsub- 
ara frequency sum and an analytic continuation to real 
frequencies need to be performed. The leading tempera- 
ture dependence of the inelastic scattering rates actually, 
comes from a branch cut in this analytic continuation.!!^! 
It turns out that to capture this effect one needs to re- 
tain terms that were neglected in the derivation of the 
integral equations in Sec. [nj Alternatively, the inelastic 
scattering rates at T = as functions of real frequency, 
or the distance in energy space from the Fermi sur face , 
can be obtained from the results of Sec. ||, see Eqs. (3.4) 
above. Scaling theory can then be used to obtain the 
results at e = and T > 0. We will follow this second 
route. 



B. Scaling Considerations 

In this section we develop a general scaling theory for 
the physical observables near the quant um critical point. 
The explicit solution given in Sec. 



z d = 2 + In g (In b) /In b 



(3.6b) 



II D is used to identify 



the results for the critical exponents. 



There are other diffusive time scales, however. For in- 
stance, the quantity H(iQ n )£l n /N-p is dimensionally a 
frequency that defines a time scale with critical exponent 



-'d 



(3.6c) 



This is the scale dimension of the phase relaxation rate, 
Eq. (3.4t). For the other two independent critical ex- 
ponents, we pick 77, which describes the wavenumber de- 
pendence of the order parameter suscep tibility , and the 
correlation length exponent v. From Eq. (2.17b), we have 
for the former 

r j = A-d-\ng{\nb)/\nb . (3.6d) 

The exponent v was determined in Ref . || with the result 

l/i/ = d-2 + lnp(ln6)/ln6 . (3.6e) 

The scale dependent exponents shown determine the 
asymptotic critical behavior, including the leading log- 
arithmic corrections to scaling. They do not include 
logarithmic terms that are less leading than the log-log- 
normal dependence due to the function g(ln&). In par- 
ticular, simple powers of logarithm s w ould correspond to 
terms of order In In 6/ ln6 in Eqs. (3.6). 



2. Thermodynamic quantities, and the density of states 



1. Critical exponents 

From the explicit solution quoted in Sec. [□] above, we 
obtain various critical exponents. As has been pointed 
out in I, there are multip le dy namical exponents. This 
i s also obvious from Eqs. ( 2.14 ) in conjunction with Eqs. 
(2.17). Furthermore, the logarithmic corrections to scal- 
ing that characterize the solution of the field theory mean 
that the asymptotic critical behavior is not given by sim- 
ple power laws. A convenient way to account for that 
is to write the critical behavior as power laws with scale 
dependent critical exponents. For instance, the critical 
time scale, which determines the dynamics of the para- 
magnon propagator, is given by a dynamical exponent 



z c — d + lng(lnfo)/ In b 



(3.6a) 



with b an arbitrary renormalization gr oup length scale 
factor. To see t his, consider Eq. (3.5) in conjunction 
with Eq. fl2.17bj ), which show that the frequency in 
the paramagnon scales like SI ~ k 2 Z) s (k, il = 0) ~ 
|k| d /g(ln(fcp/|k|)). In addition, there are diffusive time 
scales with power 2 and various logarithmic corrections. 
For instance, the frequency Q in the dressed diffuson, Eq. 
( |3.3a ), defines a time scale with a critical exponent 



We start by considering the thermodynamic properties 
near the phase transition. They can all be obtained by 
a scaling ansatz for the free energy as a function of t, T, 
and h, with h the magnetic field. Two key ideas will be 
used. First, the existen ce of two essentially different time 
scales, see Sec. Ill B 1 , implies that the free energy den- 
sity, /, should consist of two scaling parts. The second 
idea is that h also represents an energy scale (namely, the 
Zeeman energy) and it therefore scales like the frequency 
or the temperature. 

Taking all of this into account, the natural scaling 
ansatz for / is 



f(t,T, h) = b- (d+z ^ fi{tb l / v ,Tb z %hb z ") 
+b -(d+~z d ) h(tb l ' v ,Tb id ,hb z °) , 



(3.7) 



with fi and ji scaling functions. Note that a priori there 
is no reason for the temperature argument of ji to be 
given by the diffusive time scale with the dynamical ex- 
ponent Zd rather than by the critical time scale with the 
dynamical exponent z c . It requires explicit calculations 
to see that for some quantities, e.g., the magnetization, 
the diffusive temperature scale is the relevant one. In Ap- 
pendix ^| we show that the equation of state contains only 
the diffusive temperature scale, and our scaling ansatz for 
$1 reflects this feature (see also Ref. 0). In addition to 
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the T and h dependences in Eq. (3.7) with scales deter- 
mined by z c , there are subleading dependences involving 
the diffusive time scale that we suppress. The magne- 
tization m, specific heat C, and spin susceptibility Xs, 
respectively, are given by 



Fermi liquid value of N. AN can be related to a correla- 
tion function that has scale dimension — (d — 2).E_3 This 
implies the scaling form 



m = dj/dh , 
C = -Td 2 f/dT 2 
Xs =d 2 f/dh 2 . 



(3.8a) 
(3.8b) 
(3.8c) 



AN(t,e,T) = &-< d - 2 ) 

x f N (tb d - 2 g(ln b), eb d g(ln b),Tb d g{\n b)) 



(3.12a) 



Generalized homogeneity la ws fo r m, 7c = C/ T and Xs 
are obtained by using Eq. (3.7) in Eqs. (|3.8|). Substi- 



Notice that the scale dimension of AiV is minus that 
of t, modulo the logarithmic corrections to the latter. 
This leads to a resonance in the RG flow equations for 
AN, which in turn leads to an additional logarithmic 
dependence of AN on t, see Appe ndix Anticipating 



, ,• ,i , i . — i Srsi c j that logarithm, we generalize Eq. p. 12a) to 
tutmg the exponent values given by Eqs. (3.6), we find 6 ' 6 H v l -j ; 



m(t,T,h) = b- 2 f m (tb d - 2 g(lnb),Tb 2 g(lnb),hb d g(\nb)) , 

(3.9a) 

lc (t,T, h)=gQnb) 

x / 7 (tb d - 2 g (In 6) , Tb d , g (In b) , hb d g (In b) ) , 

(3.9b) 

Xs {t,T, h) = b d - 2 g{\nb) 

x f x (tb d - 2 g(\n b) , Tb 2 g(ln b) , hb d g(\n &)) , 

(3.9c) 

with f m , / 7 , and f x scaling fun ction s. With suitable 
choices of the scale factor b, Eqs. ( |3.9|) imply 



m(t,0,0) = / m (l,0,0) 



m(0,0,h) = f m (0,0,1) 



tg 



1 m 1 



d- 2 t 



2/{d-2) 



2/d 



7 c(0,T,0) = / 7 (0,l,0) 3 Qlni 

x s (t,o,o) = / x (i,o,o)r 1 , 



(3.10a) 

(3.10b) 
(3.10c) 
(3.10d) 



for the leading asymptotic c ritical b ehavior in the sense 
explained at the end of Sec. Ill B 1 . For the critical ex- 
ponents (3, 7, 8, and a, defined by m oc t 13 , \s oc t 1 , 
m oc h 1 / 6 , and 7c oc J 1 ^( 1 + Q ) (the latter is a general- 
ization to zero-temperature transitions of the usua l def - 
inition of the exponent a), we obtain from Eqs. ( 3.10| ) 



(i = 2v , 

7 = 1 , 
S = z c /2 
a = —d/z c 



(3.11a) 
(3.11b) 
(3.11c) 
(3. lid) 



with v and z c from Eqs. (3.6). 

Next we consider the single-particle density of states. 
We define AA^ — N — Np, where Np is the disordered 



AN(t, e, T) = const, x t g(ln b) In b 

+b- ( - d - 2 1f N (tb d ~ 2 g(\nb),eb d g(\nb),Tb d g(\nb)) . 

(3.12b) 

This implies 



AT(0,e,0) = N F 



CN 



■f] 



+ 0(e (d-2)/2 



1, £F 

-In — 

a e 



(d-2)/d 



(3.13a) 



and 



N(t, 0, 0) = N F 



l + d N tg 



d- 2 t 



In- 



(3.13b) 



with constants cjv and djy. It should be pointed out, 
however, that one sh ould not take the tg(\nl/t)\nt be- 
havior in Eq. (3.13b) too seriously. The reason is that 
the log- log-normal factor g(hil/t) may have multiplica- 
tive simple-log corrections that the asymptotic solution 
of th e fiel d theory is not sensitive to, see the remark after 
Eq. O. 



3. Transport coefficients 

The scaling theory for the transport coefficients, a, D, 
and D s can be presented in several different ways. Here 
we give two arguments for the scaling part of a. The 
first one starts with the fact that the conductivity is a 
charge current correlation function whose scale dimen- 
sion with respect to the quantum magnetic fixed point 
is expected to be zero. That is, a neither vanishes nor 
diverges at this quantum critical point. However, a will 
depend on the critical dynamics, since the paramagnon 
propagator enters the calculation of a in perturbation 
theory, see Appendix [b|. The critical correction to the 
bare or background conductivity further depends linearly 
on the leading irrelevant operator, which we denote by u. 
The latter is related to diffusive electron dynamics, and 
one therefore expects the scale of u to be the same as in 
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disordered Fermi liquid theory, namely, [u] — — (d — 2). El 
Scaling arguments then suggest a generalized homogene- 
ity law 

a(t,T,n) = f a (tb d - 2 g(lnb),Tb d g(\nb),nb d g(lnb), 

ub-( d -V) 

= const. + b~ {d ~ 2) f a (tb d - 2 g(lnb),Tb d g(lTib), 
nb d g (In b)) . (3.14a) 

Again, there is a resonance condition that leads to a sim- 
ple logarithm in the i-dependence of a. This is due to 
the scale dimension of u being minus that of t, modulo 
logarithmic corrections. As in the case o f the density of 
states, we therefore generalize Eq. ( [3.14a ) to 



a(t, T, fi) = const. + const, x t g(\n b) In b 

+b -(d-2) (Md-ZgQn 5) t Tb d g(\n b) , VLb d g{\n b)) . 

(3.14b) 



This scaling relation yields 



a(0,T,0) =a 



l + c CT 



T 1 e F 
7> ff U ln T 



-0(T 



{d-2)/2\ 



and 

a(t,0,0) 



1 + d a t g 



d- 2 t 



hi- 



(d-2)/d 

(3.15a) 

0(t) , 
(3.15b) 



with c a and d a constants. The same caveat as given after 
Eq. ( pTT3h| ) applies. 



An alternative argument that gives Eqs. (3.15) is to 
assume that a consists of a background part that does not 
scale, and a singular part, 8a, that does. In fundamental 
units, [Scr] = — (d — 2). This suggests the scaling form 



6a{t,T,Q) =b- (d - 2) 

xf a (tb d - 2 g(\nb),Tb d g(lnb),nb d g(\nb)) 



(3.16a) 



Here the dependence on u is already implicitly taken into 
account, so we dropped the explicit dependence. Taking 
into account the resonance between the scale dimension 
of 5a and t, see Appendix |d|, gives 



Sa(t, T, il) = const, x t g(ln b) In b 

+b- {d - 2) f a {tb d - 2 g{\nb),Tb d g{\nb),nb d g{\nb)) 



(3.16b) 



This also yields Eqs. (|1|). 

The diffusion coefficients D and D s have dimensions of 
length squared divided by time. Since there are two time 
scales, there are two possible scale dimensions for the dif- 
fusion coefficients. is the diffusion coefficient for the 



order parameter fluctuations, so one expects the critical 
time scale to apply, while D is the quasi-particle diffu- 
sion coefficient, so the diffusive time scale is appropriate. 
This leads to homogeneity laws 

D s (t,T,n) = [&-( rf - 2 Vs(ln6) 

xfD s (tb d ~ 2 g(\nb),Tb d g(lnb),nb d g(\nb)) , (3.17a) 

D(t,T,Q) = [gilnb)}- 1 

xf D (tb d - 2 g{\nb),Tb d f{\nb),nb d g(lnb)) . (3.17b) 

These two results are consistent with the fact that the 
conductivity is noncritical to leading order, and scales 
like a ~ D s x s ~ DH. Indeed, the above results for a, 
D s , and D can be used to obtain scaling results for Xs 
and 7c oc H , which justify our free energy considerations 
in Sec. Ill B 2 above. 



4- Relaxation rates 

The phase breaking rate r^ 1 and the quasiparticle de- 
cay rate Tgp are given by Eqs. (3.4). By comparing 
Eq. (|3~4b|) and Eq. (|3~3a| ), we see that r pt scales like 
a wavenumber squared (recall that G is not singularly 
renormalized and hence does not scale), and therefore 
has a scale dimension [t~^] = 2 with no logarithmic cor- 
rections. This observation leads to the homogeneity law 

r; h \t,e,T) = b- 2 

x / ph (^ d - 2 5 (ln b) , eb d g(\n b) , Tb d g{\n b) , ub~ {d - 2) ) . 

(3.18a) 

The leading irrelevant variable u represents interaction 
effects that are necessary for any dephasing. The rate is 
therefore linear in u, and we can write 

r^{t,e,T)=b- d 

xf ph (tb d ' 2 g(\nb),eb d g(lnb),Tb d g(\nb)) . (3.18b) 

At criticality, we find 

7^(0, e,0) =c ph e 9 Q In ^ +0{e d ' 2 ) , (3.18c) 

with c p h a constant. 

The quasi-par ticle relaxation rate is given by the ratio 
H"/H', see Eq. (3.4a). The scaling properties of Tgp thus 
follow from those of H, or jc, Eq. (3.9b). Explicitly we 
find 



t qp(* = Q,e,T = 0) = c Q p eln 



hi 



/In 



(3.19) 



10 



with cqp a constant. For e <C cf, £ 7^ we have asymp- 
totically 



'QP 



(3.20) 



To obtain Eq. ( 3.2C ) we have used the well-known fact 
that at the disordered Fermi-liquid fixed point the re- 
laxation rates are proportional to e d / 2 . In terms of our 
scaling arguments this result is rederived in Appendix ^. 



IV. DISCUSSION, AND EXPERIMENTAL 
RELEVANCE 

This paper completes our discussion of various aspects 
of the ferromagnetic phase transition in low-temperature 
disordered itinerant electron systems. We have given the 
exact solution for the magnetic critical behavior near the 
quantum phase transition from a paramagnetic metal to 
a ferromagnetic metal. We have also determined the crit- 
ical behavior of a number of other relevant physical vari- 
ables near the transition, in particular the electrical con- 
ductivity and the tunneling density of states. In addition, 
we have established several connections between previ- 
ously formulated theories. We conclude with a discussion 
of general aspects of our results, and their experimental 
relevance, starting with the former. 



A. General discussion 

One important aspect of the present analysis is the 
establishment of connections between various theoreti- 
cal formulations of the quantum ferromagnetic transition 
problem. First, we related our previous nonlocal order 
parameter field theorycl to the local coupled field the- 
ory involving both the fermion density fluctuations and 
the magnetic order parameter fluctuations that was for- 
mulated in I. The two theories yield the same critical 
behavior apart from logarithmic corrections to power- 
law scaling that were missed in Ref. ||, and are correctly 
given by the coupled local theory. Second, we have un- 
ambiguously related earlier work on the disordered in- 
teracting electron problem, involving runaway renormal- 
ization group flow, to the quantum ferromagnetic phase 
transition.E3 This connection was actually made earlier 
in Ref. ||. However, as explained in the Introduction, our 
argument then was not complete. Indeed, even though 
Ref. U correctly obtained the critical behavior, it failed 
to identify the nature of the phase transition as the fer- 
romagnetic one. In hindsight, this is surprising, given 
that scaling theory was used to correctly identify the or- 
der parameter exponent /? = 2v. The problem was that 
Ref. H was formulated solely in terms of fermionic num- 
ber and spin density fluctuations at the Fermi surface, 
so that the behavior of quantities that involve electrons 
far from the Fermi surface, such as the magnetization, 
was not obvious. Indeed, we argued that although the 



above exponent equality was formally valid, the magne- 
tization was actually zero in the ordered phase since the 
scaling function had a zero prefactorO This argument 
was incorrect because it failed to take into account that 
electrons away from the Fermi surface are also ordering. 
This, in turn, leads to a nonzero scaling function. 

The difficulties interpreting the theory put forward 
in Ref. |^ notwithstanding, it is very remarkable that 
this nonlinear sigma model formulation of the problem 
yielded the correct result, since it was not geared at all 
towards describing ferromagnctism. The focus on degrees 
of freedom near the Fermi surface mentioned above is one 
reason, and another one is the fact that the sigma model 
is derived by expanding about the paramagnetic metal 
fixed point, so it is not obvious why it is capable of de- 
scribing a critical fixed point. rlhis is actually a general 
question about sigma models, E3 and the answer is only 
incompletely known. Indeed, periodically even the capa- 
bility of the O(N) nonlinear sigma model to qualitatively 
correctly describe_the Heisenberg transition in d = 3 has 
been questioned. Ej 

The connection between the runaway renormalization 
group flow encountered in low orders of a loop expan- 
sion and ferromagnetism is particularly interesting in 
two-dimensional systems because of recent experiments 
that show either metallic or metallic-like behavior in Si 
MOSFETs and other materials,E-3 and even more recent 
ones that show that this behavior happens near a-quan- 
tum phase transition to a ferromagnetic stateO The 
connection between ferromagnetism and two-dimensional 
metallic-like behavior is not obvious, but the observa- 
tion of a ferromagnetic phase in d = 2 is consistent 
with our identification of the runaway flow behavior with 
ferromagnetism. The .same connection was more re- 
cently made by others.t^l It is also interesting to note 
that a nearby ferromagnetic phase in disordered sys- 
tems is favorable tp-jan exotic type of even-parity, triplet 
superconductivity.EZI After the experimental observation 
of 2-d metallic-behavior, this was proposed as a possi- 
ble explantionE3 Proposals of superconducting or other- 
wise exotic phases as the explanation for the observations 
are bolstered by the conclusion that conventional metal- 
lic behavior in a two-dimensional ferromagnetic system 
is unlikely to occur.Ej Even if the observed metallic or 
pseudo-metallic phase is unrelated to superconductivity, 
the presence of a ferromagnetic phase makes the existence 
of a triplet superconducting phase nearby more likely. 



B. Experimental consequences 

Most of the results of the present paper can be directly 
checked by experiments, at least in principle. For in- 
stance, the pressure tuned ferromagnetic [transitions ob=. 
served at very low temperatures in MnSO and UGe2E2l 
provide examples of systems where the quantum critical 
point is directly accessible. These experiments were done 
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on very clean samples, where the ferromagnetic transi- 
tion at low temperatures is ,©f first order, in agreement 
with theoretical predictions.^ However, upon introduc- 
ing quenched disordeppne expects the transition to be- 
come of second order and the current theory to apply. 

The critical behavior predicted for the thermodynamic 
quantities is markedly different from_the mean-field ex- 
ponents predicted by Hertz's theory.El For instance, the 
predicted value of the magnetization exponent in d = 3, 
(3 = 2 with logarithmic corrections, is very different from 
both the mean-field value /?mf = 0.5 and the 3-d classi- 
cal Heisenberg value /3h ~ 0.37. One important remark 
in this context is that the logarithmic corrections will, 
over any realistically achievable range of t- values, mimic 
a power, so that the observed value of /3 should be ex- 
pected to be smaller than 2. Similarly, the correlation 
length exponent v will be equal to 1 with logarithmic 
corrections in d — 3. However, these exponents may 
be hard to measure directly, especially near a quantum 
phase transition that must be triggered by a non-thermal 
control parameter that is more difficult to accurately vary 
than the temperature. It is therefore important that the 
values of the critical exponents are also reflected in the 
behavior of the tunneling density of states, and the elec- 
trical conductivity, across the transition. Even though 
these observables do not show any leading critic al be - 
havi or, th e leading corrections expressed in Eqs. (3.13) 
and ( 3.15 ) reflect the values of v and z, and they should 
be easier to measure than the critical behavior of, say, the 
magnetization. For instance, our prediction for the tun- 
neling density of states in d — 3 is as follows. Far from 
the transition, it will show the well-known square-root 
anomaly as a function of the .bias voltage that is charac- 
teristic of disordered metalsJia Near the transition, the 
voltage region that shows the square-root behavior will 
shrink, and outside of it a region of cube-root behavior 
will appear, until at criticality the b ehavior is given by 
the e 1 / 3 behavior shown in Eq. ( 3.13a ). The same discus- 
sion applies to the co nductivity as a function of tempera- 
ture, see Eq. ( [3.15a ). Again, the logarithmic corrections 
to scaling will manifest themselves in a real experiment 
as an effective power smaller than 1/3. 

In this context we also come back to the relative values 
of the coefficients and ad-2 in Eq. ( |2.3a| ). The leading 
critical behavior is due to the |k| d ~ 2 term whose coeffi- 
cient is a,d~2, as can be se en, f or instance, from Eq. ( p32| ) 
in conju nction with Eq. (2.6c). Since ad-2 — 0(1 /kpt), 
see Eq. ( [2.3a| ), this implies that the leading effects will 
be manifest only for sufficiently strong disorder, or, if 
we scale the wavenumber with the correlation length £, 
for sufficiently large £ at fixed mean-free path £. Since 
ad = 0(1), the nonanalytic term will dominate for £ > £ 
or, using v = 1, for t < l/kpt. Typical values of the dis- 
order results in mean-free paths £ ~ 10/kp. For such a 
value, the our leading results will apply everywhere in the 
critical region. For less disordered samples, their region 
of validity will be correspondingly narrower. 



C. Conclusion 

In conclusion, we now have a complete theory for the 
quantum critical behavior of disordered itinerant ferro- 
magnets in d > 2, including the exact values of the 
critical exponents, the leading logarithmic corrections to 
power-law scaling, and the relations between various the- 
oretical approaches to the problem. Specific predictions 
for the behavior of all important observables allow for 
a direct experimental test of this theory. However, in 
d = 2 there remains a puzzling discrepancy between ex- 
isting theory and observations. The latest experimental 
evidence is for a transition, with increasing electron den- 
sity, fr-epi a paramagnetic insulator to a ferromagnetic 
metaljEa while there is no theory that can account for 
a metallic state, ferromagnetic or otherwise, in d = 2. 
In particular, it has recently been shown that ferromag- 
netic fluctuations in d = 2 do not produce, a metallic 
state within a perturbative RG treatment fB ruling out 
a possible mechanism for a metal-insulator transition in 
d = 2. This state of affairs has recently been reviewed in 
Ref. M 
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APPENDIX A: ADDITIONAL 3-POINT AND 
4-POINT VERTICES 

Starting at one-loop order, the RG generates vertices 
that are not in the effective action. An example is the 
3-point vertex shown in Fig. ^. In contrast to the vertex 




FIG. 9. A one-loop diagram that generates a new 
three-point vertex. Parts of the diagrams shown in Fig. ^ 
also contribute to this vertex. 

C2, all external legs in this diagram carry the same replica 
index. To see the physical meaning of this term, we in- 
tegrate out the 6-field to arrive at an effective 4-point q- 
vertex that is completely diagonal in replica space. This 
corresponds to a four-body interaction term which in 
imaginary time space must have the form 



dr (u s (t)) £ 



(Ala) 
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with n s (r) an electron spin density in imaginary time rep- 
resentation. Performing a Fourier transform, and making 
use of the isomorphism between density operators and q- 
matrices that has been explained in I (see also Ref. 15), 
this corresponds, in schematic notation, to 



T 3 | dxJ^Q^x)) 4 



(Alb) 



Here Q(x) is the matrix field from I Eq. (2.8), and we 
have suppressed frequency labels and frequency sums for 
clarity. A Hubbard-Stratonovich transformation to rein- 
troduce the 6-field then leads to a bq 2 vertex of the struc- 
ture 



c 2 T 3 / 2 



dx^V Q (x) 



? Qa (x) q QQ (x) 



(A2) 



with C2 a coupling constant. This vertex thus carries a 
higher power of the temperature than the one with cou- 
pling constant C2. By using Eqs. ( f2.ll ) to estimate the 
behavior of the diagram shown in Fig. ^| we find that 
C2 diverges for d < 4 in the long-wavelength and small 

with A the infrared mo- 
The scale dimension of 



frequency limit, ci — c 2 Af 
mentum cutoff from Sec 



[IB2 



C2 is therefore smaller than that of C2, 
[B 2 ]<[c 2 ]-(d-2) . 



(A3) 



This term is therefore irrelevant for the critical behavior. 

Similarly, the RG generates four-point vertices that are 
not in the effective action. For instance, diagram (a) in 
Fig. |l(] regenerates the two-body interaction that was 
shown in I to be irrelevant and therefore dropped. Dia- 
gram (b) is a four-body interaction of the same type as 
discussed above in connection with Fig. [| All of these 
terms, and similar ones not shown in Fig. |l0|, are thus 
RG irrelevant and can be safely neglected. 



(a) 



(b) 



FIG. 10. One-loop diagrams that generate (a) a two-body 
and (b) a four-body interaction. 



Using Eqs. (2.7), d2_J), this yields 



Q_ 1 

AN(e F + e) = - G 2 NpK t — £ T £ (V, 



m<0 



xK-m(p) 



(B2) 



In Sec. H we have seen that self-consistent one-loop the- 
ory for the propagators T> and M. is exact. Using the 
resulting dressed propagato rs in Eq. (B2) and doing the 
integrals, one obtains Eqs. ( 3 . 1 3| ) . 

Similarly, if we use the perturbative result for the con- 
ductivity given by Eqs. (3.6) of I and the exact propa- 
gators derived in Sec. || of the present paper, then Eqs. 
( |3~15l ) for a are obtained. 

These two results buttress the scaling arguments given 
in Sec. 



Ill 



Note that these perturbative calculations 
are equivalent to taking the standard Altshuler-Aronov 
perturbative resultsJIS and using propagators that are 
appropriate near the magnet ic q uantum phase transition. 

The classical limit of Eq. ( |B2| ) for the density of states, 
or the equivalent result for the conductivity a, can be re- 
lated to the established result for the critical behavior 
of the conductivity at a Heisenberg critical point £3 To 
see this, note that at fini te temperatures the diffusion 
propagator, T> in Eq. (B2), has a mass due to inelastic 
scattering processes. For the leading critical behavior, T> 
can therefore be replaced by a constant. Also, in the clas- 
sical limit, the frequency sum in Eq. (B2) turns into an 
integral over all frequencies, since in this limit the Boltz- 
mann weight that restricts the frequency sum is absent. 
The net result is that this contribution to AN or Ac 
is proportional to a magnetization-magnetization corre- 
lation function that is local in space and time, 



AN (x Act oc (M 2 (x,r)) 



(B3) 



The correlation function on the right hand side of Eq. 
( B3 ) is essentially the magnetic energy density, and hence 
scales as t l ~ a , with t the distance in temperature space 
to the classical phase transition, and a the usual specific 
heat critical exponent. This gives the result of Ref. p2l 



AN ~ Act ~ t 1 



(B4) 



APPENDIX B: PERTURBATIVE RESULTS FOR 
THE DENSITY OF STATES AND THE 
CONDUCTIVITY 

In this Appendix we show that the scaling results ob- 
tained in Sec. Ill are consistent with perturbation theory 
for N and a. 



From Eq. (3.1) we find the one-loop result for AN as, 
AN(e F + e) 

z ' — ■ ' — • — 

j3 m i,r 

XrCm( X )>iu%^e+iO . (Bl) 



APPENDIX C: EQUATION OF STATE 

In this Appendix we show how to generalize the self- 
consistent equations given in Sec. || for correlation func- 
tions in the paramagnetic phase to the ferromagnetic 
phase. In this way we will derive and validate the scaling 
argument given in Sec. Ill for the magnetization m as a 
function of t. 

To simplify the discussion we first ignore the logarith- 
mic corrections to scaling. Then, in the paramagnetic 
phase, the diffusion coefficient D is simply a number, 
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and the equation for D s is given by Eq. (2.14a) with 
D(u>) replaced by that number. The propagators in the 
ordered phase have been derived in Ref. |lj. Assuming 
that the magnetization is in the z-direction, we first need 
to specify whether longitudinal or transverse spin den- 
sity correlations will be considered. In the ferromagnetic 
phase, the transverse spin-density fluctuations become 
propagating Goldstone modes or spin wave excitations, 
while the longitudinal ones remain diffusive. For the lon- 
gitudinal spin-density mode, the left hand side of Eq. 
( p. 14a ) therefore still describes a diffusion coefficient. It 
is easily shown that in this case the diffusion poles on 
the right hand side of this equation are cut off by the 
magnetization. The frequency corresponding to this cut- 
off is a "cyclotron" frequency that scales like uj c ~ \m\. 
The explicit generalization of these propagators to the 
ferromagnetic phase is given by Eqs. (3.10) o f Ref. |§ 
Using these results, the generalization of Eq. (2.14a) to 
the ferromagnetic phase is, at zero external frequency 
and wavenumber, 



iG_ 
2V 



E 



du> ■ 



(p 2 — iuj/D + c|to|) 5 



(CI) 



with c a constant. Carrying out the integrals yields, for 
2 < d < 4, 



-ci\t\ 



C2 TO 



(d-2)/2 



(C2) 



Here \t\ = — t > is the distance from criticality in the 
ferro mag netic phase, and c\ and C2 are constants. Equa- 
tion (C2) implies that to scales as 



TO - |i| 2 /(^ 2 ) 



(C3) 



Equation (C3) is consiste nt with the scaling result for the 
critical exponent /3, Eq. ( 3.11a ), apart from logarithmic 
terms. 

The logarithmic corrections to the exponent (3 can be 
understood as follows. As noted in Sec. [II Bl, the cor- 



relation length exponent v has leading logarithmic cor- 
rections, while the critical exponent 7 does not.tl Within 
the integral equation approach, and in the paramagnetic 
phase, this manifests itself in the following structu re o f 
the renormalized paramagnon propagator (see Eqs. (3.5) 
and (p7b|)), 



t + \k\d-y g (\n(k F /\k\)) + \n n \/v 



(C4) 



That is, the term proportional to |k| d ~ 2 carries leading 
logarithmic corrections, while the term t does not. In 
Eq. (C4) we have left out all constants for clarity. If 
one scales the wavevector with the correlation length, 
then this structure produces the logarithmic corrections 
to the exponent v. In the ferromagnetic phase, the |k| d ~ 2 



nonanalyticity is cut off by a magnetic length or cyclotron 
radius, £ m oc 1/m 1 / 2 . This means that the |k| d_2 gets 
replaced by m^ d ~ 2 ^ 2 , see Eq. (|C2|). The net result is 
that, again, t has no leading logarithmic correction, while 
TO,(d-2)/2 does have one. Sc aling to with the appropriate 
power of t then yields Eq. ( 3.11a ) for (3. 



APPENDIX D: LOGARITHMIC CORRECTIONS 
TO SCALING 

Wegner0 has given a classification of logarithmic cor- 
rections to scaling. The first class consists of simple log- 
arithms that arise due to resonance conditions between 
scale dimensions. In the present context, such a reso- 
nance occurs between the scale dimensions of the leading 
correction to the single-particle density of states and the 
relevant variable t. For clarity, let us neglect the more 
complicated logarithms that are embodied in the function 
g(\nb) for the time being. The one-loop flow equations 
for these two quantities then are 



dAN 
dlnb 
dt 
dhd> 



(d- 2)AiV + const, x t 
(d-2)t . 



(Dla) 
(Dlb) 



The general solution of the homogeneous equation for 
AN is 



(AN) hom (b) = (AN)(b = l)b 



d-2 



(D2) 



This has the same ^-dependence as the inhomogeneity, 
t(b). Consequently, the solution of the inhomogeneous 
equation is 



(AN)(b) = [(AN)(b = 1) + const, x t(b ■ 



l)ln&] b d - 2 . 

(D3) 



The resulting logarithm has been taken into account 
in Eq. (3.12b). (Notice that the physical quantity is 
AN(b = 1).) The same mechanism is at work for the 



conductivity, and this is reflected in Eq. (3.14b). 

Wegner's second mechanism is due to marginal opera- 
tors, and it can lead to arbitrary functions of logarithms. 
In our case, C2 acts as an effectively marginal operator, 
as has been explained in I, and this leads to the log-log- 
normal factors we denote by g(ln6). 



APPENDIX E: RELAXATION RATES IN A 
FERMI LIQUID 

Here we illustrate how to obtain Schmid's r esultS for 
the relaxation rates in a disorder ed F ermi liquid from the 
scaling theory developed in Sec. III. 

At a disordered Fermi liquid fixed point, the dynamical 
exponent is z = 2, reflecting the diffusive dynamics of the 
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quasiparticles, and the leading irrelevant variable, which, 
we denote by it, has a scale dimension [u] — — (d — 2).t£l 
Since ~ Tgp both are dimensionally frequencies or 
energies, they scale the same way and we have for either 
rate a homogeneity law 



T -i( e , T) = b- 2 f T (eb 2 , Tb 2 , ub~^) 



(El) 



The dependence on u arises from the electron-electron 
interaction terms that lead to r _1 ^ in the first 
place, and therefore the scaling function has the prop- 
erty / r -i (1, 0, x) oc x. The explicit dependence on u can 
therefore be eliminated by writing, instead of Eq. (El), 



T-'(e,T) = b- a Meb 2 ,Tb 2 ) 
In particular, we have@ 

r- 1 (e,0) = / T (l,0)e d / 2 
which we used to derive Eq. ( |3.2C| ). 
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(E2b) 
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